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SUMMARY 

A method of obtaining the flow of a nonviscous compressible fluid 
past arbitrary compressor or turbine blades between two neighboring sur- 
faces of revolution is presented. The equations of continuity and 
motion obtained for such flow are combined into a nonlinear second- 
order differential equation in terms of a stream function defined for 
such flow. Numerical solution by the use of differentiation coeffi- 
cients for unequally spaced grid points is suggested. Means for satis- 
fying the boundary conditions outside the channel and solutions by the 
relaxation method with manual computation and by a matrix method with a 
large-scale digital computing machine are described. 

Satisfactory results were obtained by the method described in the 
investigation of the detailed flow variation of a compressible fluid 
past typical hi^-solidity, highly cambered thick tvirbine blades on a 
cylindrical surface. The variations in fluid properties across the 
channel appeared to be representable by a second-degree f\mctlon. The 
mean streamline approximately followed the shape of the mean channel 
line of the cascade and had lower curvature. The variation of specific 
mass flow along the mean streamline followed the trend of the varia- 
tion in the channel width. In general, the variation in specific mass 
flow was significantly higher than that given by the ratio in channel 
width and the effect extended outside the channel. The velocity distri- 
bution around the blade obtained in the theoretical calculation compares 
very well with experimental values. 


INTRODUCTION 

A basic aerodynamic problem of tiorbojet and turbine-propeller 
engines is the flow of compressible fluid past a series of blades in 
circular sirrangement . In axial -flow type turbomachines, if the blades 
are relatively short in the radial dimension and are bounded by cylin- 
drical walls, the theoretical flow passing through the blades is usually 
computed on the basis of two-dimensional flow on a cylindrical surface. 
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which is then developed into a plane for convenience of calculation. A 
number of methods have been proposed to obtain the theoretical flow 
through such a given cascade of airfoils, including analytical methods 
using conformal mapping, interference technique, and Prandtl-Glauert or 
Chaplygin-K^mln-Tsien approximation for compressible flow (references 1 
to 6), graphical procedure (references 7 to 10), and other mechanical 
and electrical devices (reference lO) . For radial- and mixed-flow type 
turbomachines, the available methods consider inconpressible flow in a 
radial plane (references 7 and 11 to 16) or compressible flow on a coni- 
cal surface (reference 17). However, in current axial-flow turbo- 
machines, the hub and casing walls may be either tapered or ciurved, 
which causes the fluid to flow on a noncylindrical surface; and in cen- 
trifugal machines, the flow surfaces are usually quite curved toward 
the inlet (see fig. l) . A method of analysis is therefore developed at 
the NACA Lewis laboratory to analyze the two-dimensional compressible 
flow for the fluid between two neighboring surfaces of revolution in 
these turbomachines (fig. l) . The equations of continuity and motion 
for irrotational absolute flow are first obtained for such flows and are 
then combined into a nonlinear second-order partial differential equa- 
tion. Because the change of fluid properties passing through turbo- 
machine blades is, in general, large and the shapes of surface and 
blades are arbitrary, nimierical solution by the finite-difference 
approach is suggested. 

Solving the flow throu^ a cascade of blades of arbitrary camber 
and thickness by the finite -difference method involves two main diffi- 
culties; (l) the curved boundaries formed by the blade surface, and 
( 2 ) the large number of grid points necessary to cover the whole flow 
region. The first difficulty is removed by the recently available' dif- 
ferentiation formula and coefficients for unequally spaced grid points 
(reference is) . With these formulas and coefficients, the curved bomid- 
ary can be handled in the same manner as a strai^t boundary. The 
second difficulty can be reduced by the use of higher -degree polynomial 
representations (references 18 and 19) . Furthermore, if a modem high- 
speed large-scale digital machine is available, the set of simultaneous 
difference equations in the stream function can be very quickly solved 
by the matrix process given in reference 18, and consequently, suc- 
cessively improved solutions for compressible flow can be obtained in a 
reasonably short time. Without such a machine, it is desirable to 
obtain an approximate solution for compressible flow by either graphic, 
mechanical, or other approximate methods and to use the relaxation 
method (references 20 and 21) for the final improvements. 

In t\arbomachines with a relatively large radial dimension, curved 
hub and casing walls, or if designed on velocity diagrams other than the 
free-vortex type, the radial flow assumes primary importance. For such 
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turboma chines^ methods of flow analysis have been proposed that take 
into account the radial flow but consider only the variation of the 
fluid properties on a mean stream surface passing between two blades 
(reference 22) . When applying such methods of anailysis to turbomachines 
with thick blades, it is necessary to have some knowledge of the effect 
of area reduction due to blade thickness on the mean flow and the rela- 
tion between the blade shape and the shape of the mean streamline at 
various radii. For this reason, the method developed is applied to 
investigate the detailed compressible flow in a typical turbine cascade 
with highly cambered thick blades. The theoretical velocity on the 
blade is compared with experimental data obtained at the NACA Lewis 
laboratory. 


SYMBOLS 

The following symbols are used in this report: 
a,b,c points on streamline 
B differentiation coefficients 

C constant 

E error 

H total enthalpy based on absolute velocity 

h static enthalpy 

I h + I W2 - I o)2r2 

L blade length projected on z-axis 

I orthogonal coordinates on mean surface of revolution, (fig. l) 

M total mass flow passing throiigh space between two neighboring 

blades 

n degree of polynomial 

P pitch or spacing 

p static pressure 

r radial distance from the axis of machine^ (fig* l) 
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t 

W 

y 

z 

a 


blade thickness in pitch direction 

velocity relative to blade 

rep 

distance along axis of machine, (fig. l) 


flow angle in 2cp-plane, 


tan"l 


1 Wep 
r ^ 


y ratio of specific heats 

S grid spacing 

p mass density 

a tan-J- (fig. 1) 

T normal thickness of stream filament of revolution, (fig. l) 

^ stream function 

o) angular velocity of hlade 

Subscripts: 

i inlet 


j,k grid points where differentiation coefficients are applied 

I,cp meridional and circumferential component, (fig. l) 

m mean streamline 

n degree of polynomial 

p pressure surface 

s suction surface 

T total, or stagnation, state 
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Superscripts: 

grid points 

1^2 indicate first and second derivatives, respectively 


METHOD OF SOLUTION 

The steady isentropic flow of a nonviscous compressible fluid along 
a stream filament of revolution described by the orthogonal coordinates 
I and cp on the mean surface of revolution and a varying normal thick- 
ness of the filament T (fig- l) is governed by the following equa- 
tions of continuity and irrotational absolute flow and the isentropic 
pressure-density relation (see appendix for derivation of the first two 
equations) : 


S(TPWir) ^(TpWcp) ^ 

■ Tcp " = ° 


( 1 ) 


hVI. 


cp 


1 d 


r ^cp 


+ 2o) sin a = 0 


(2) 


p = C pT^ 

From equation (l) , a stream function ijr can be defined by the 
following relations: 

TPW^ - - I 

Substituting equations (4) in equation (2) results in 

+ ( sin b SlnX^Silf, 1 

^^2 V r ■ ~Sr" JWl ^2^ 


(3) 


(4) 


Slnp^t, iSln p3\lf,„ 

~W~ ^ ~5'cp'"" + 2 o)Tp sin q j = 0 


( 5 ) 
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In the n-umerical solution of the problem, iteration is necessary 
for compressible flow. A network of grid points is chosen to cover the 
whole flow region, and finite-difference equations are obtained for the 
first three terms of equation (s) . The last term is taken as constant 
during the solution of the \|r values and recomputed from the improved 
values of ijr for the next solution (references 20 and 2l) . In this 
way, the numerical coefficients at each grid point are determined by 
the geometry of the problem alone. 

The value of density is most conveniently expressed throughout the 
flow region in terms of velocity and enthalpy by the use of equation (3) 
and the relation for a perfect gas: 


1 



where I is constant throughout the flow region for adiabatic, abso- 
lutely irrotational, and steady relative flow (reference 22). (For 
application to stator blades, o) becomes zero, W becomes the abso- 
lute velocity, and I becomes H. ) If the components of W as 
expressed by equations (4) ajre substituted in equation (6) , the follow- 
ing equation is 


P 

PT,i 


I + ^ 


Hi 


1 



2 Hi 


1 

i>2 \5y/ 


r-1 


( 7 ) 


In order to facilitate the evalxxation of the density from the 
\lf-derivatives, equation (?) can be rewritten in the following form: 


2 

2 = (1 - 


(7a) 


where 
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(7b) 
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and 


c> = 






(7c) 


With the relation between p and ijr -derivatives in this form, either a 
table may be constructed by computing 2 for equally spaced values of $ 
through iteration on a punch card machine, or a graph may be constructed 
by computing $ for a number of Z. Because I, oo, r, 

and T are all known values, with the aid of the table or graph, the 
value of density at any point can be easily obtained after the first- 
order derivatives of ilr are computed by using the differential coeffi- 
cients for equally and unequally spaced grid points given in refer- 
ences 23 and 18, respectively. 


The first-order derivatives of density are then computed and com- 
bined with the first-order derivatives of If the last term in equa- 

tion (5) is denoted by J and the differentiation coefficients, which 
mixltiply the \|r value at point J to give the second and first 
derivatives at point i using an n't^^-degree polynomial, are denoted by 

and respectively, the finite -difference form of eqiiation (5) 

at the grid point i can be written as 



sin a 


In T 

Sz 


^B^ 
n J 


n 


k=0 


+ ji = 0 


( 8 ) 


where and denote ^ values on the surface along I and cp 

coordinates^ respectively (figs. 1 and 2). If second-degree polynomial 
representation is chosen (n=2), three points are involved in the 
numerical differentiation in each direction^ and the center -point for- 
mula can always he used for all interior grid points. If fourth-degree 
polynomial representation is chosen (n=4)^ five points are involved 
in each direction and the off -center -point formula has to he used at 
the point next to the houndary. Two typical cases where odd spacings 
are used in the calculation of the coefficients because of a curved 
boundary are shown in figure 2. In the case where \|r^ is at j=2^ 
k=2^ the center-point formula is used in both directions, with equal- 
space I -direction coefficients obtainable from reference 23 and 
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vith cp-direction coefficients for the odd space at the boundary obtain- 
able from reference 18. In the other case where is at j=2, 

k=3^ again the equal-space center-point coefficients are used in the 
7 -direction; but in the qp -direction an off -center -point formula for 
coefficients must be used, including an odd space at the boundary, 
again obtainable from reference 18. It is desirable, where possible, to 
have the spacing near the boundary somewhat shorter than the regular 
equal spacing so that accuracy at the point next to the boundary is 
comparable with that at other points. By using a five-point system^ 
many less grid points are required in comparison with the three-point 
system for the same accuracy. (See reference 18.) In the present 
problem, the variation of is only rapid near the leading and trail- 
ing edges of the blade and the grid pattern should be determined for 
these regions first. 

After the degree of polynomial representation and the size of grid 
have been chosen, a finite -difference equation (8) can be easily written 
for each interior grid point using the differentiation coefficients for 
equal and .’nequal grid spacings. The constant value on the suction 
surface can be arbitrarily chosen and the constant \lr value on the 
pressure surface is determined by the mass flow passing through the 
channel or at the inlet as follows: 


^Cpp 




% 








TpW^r dCp = M 


(9) 


Outside the channel two reference lines, a pitch angle apart, can be 
drawn either parallel to the I -axis or parallel to the inlet and exit 
angles in the Icp-plane. (Parallelism with the inlet and exit angles, 
as shown in fig. 3, is a little more convenient for the purpose of 
drawing streamlines after the solution is obtained. ) Because the 
fluid state repeats itself for each pitch angle along the cp-direction, 
only the grid points lying between the reference lines need to be 
included in the calcvilation when the number of grid spaces per pitch 
angle is an integer, and the central five-point formula can always be 
used at every grid point between these lines. For example, at grid 
point c (fig. 3), the \lr values at a and b are obtained in terms 
of the values at i and j, which are a pitch angle apart from 
a and b, respectively, by the following relation: 

;lra = i|fi + - ,|fg) 

= 1 9 + (tp - tg) 


(10) 
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Inlet and Exit Boundary Conditions 


The solution of the present problem in terms of \lr has an inter- 
esting feature in the boundairy conditions. As noted in the previous 
section, the \lc value is constant along either the suction or pressure 
surfaces, which are fixed boundaries^ and outside the channel, there is 
no longer a fixed boundary , but there is the condition that the flow 
repeats itself for every pitch angle in the cp -direction. Sufficiently 
far upstream of the blade, an inlet station can be chosen so the bound- 
ary condition there requires the fluid state to be uniform in the 
cp-direction and prescribes the inlet angle The following two 

methods have been devised to account for this boundary condition: 

Fixed-angle method. - A part of the grid system in the Icp-plane 
near the inlet station i-i is shown in figure 4(a). In order to write 
the finite -difference equation at point a, the value of t a-'t 
point c to the left of station i-i can be obtained by using the given 
inlet angle Pi and the linear variation of t at station i-i 


The coefficients of and in equation (H) after having been 

multiplied by are added to the regular coefficients at points a 

and b; otherwise the points a and b are treated in the same manner 
as other interior points. Whether the inlet station i-i is chosen suf- 
ficiently far away from the blade is indicated by the linearity of the 
\|c variation at that station as obtained in the solution. 

StreanO i‘ne-ad,justment method . - If the first station i-i is chosen 
piif f 1 r»i enti y far from the blades, the variation of stream function to 
the left of the station i-i is linear in the circumferential direction. 
The value of the stream function, however, depends on the inlet angle. 

If solutions for a range of inlet angles are desired, they can be 
obtained by specifying a number of linearly varying stream functions to 
the left of station 1-i as fixed bo^lndary values. The streamline 
obtained in the solution then gives the value of the inlet angle. If, 
however, the solution for a certain inlet angle is desired, the stream- 
line obtained in the solution must be adjusted according to that inlet 




( 11 ) 
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angle, (for example as ab in fig. 4("b) is adjusted to position ac), 
thereby obtaining an improved set of boundary values of the stream 
functions to be used in the next calculation. This method is of course 
not as accurate and convenient as the fixed-angle method for obtaining 
a solution for a given Inl et angle but is desirable in the matrix solu- 
tion (to be discussed in the following section) because the inlet angle 
is not involved in the matrix factorization, making the same matrix 
factors usable for a range of inlet angles and Mach numbers. 

At the exit station far downstream of the blade, the same methods 
can be applied. For a blade having a sharp trailing edge, the Kutta- 
Joukowski condition can be used and the correct exit angle far down- 
stream is the one that gives the flow at the trailing edge satisfying 
that condition. For round trailing edges, either the position of the 
stagnation point can be assumed or some available empirical rules for 
the exit angle used. If the calciilation is made to compare with cer- 
tain experimental results, the measured exit angle can be used. 


Solution of Finite-Difference Equations 

With the grid system and the degree of polynomial representation 
chosen and the boimdary conditions taken into account, the problem 
remaining is the solution of the set of N linear algebraic equa- 
tions (8) written for N interior grid points. For a small number of 
solutions with a given blade, the best method is the relaxation method 
(reference 20) . A modification of this method involving the use of 
higher-order differences is suggested by Fox (reference 18). Formulas 
and the table of coefficients obtained in reference 18 enable the 
direct use of higher-degree polynomials for problems with curved bound- 
aries. For the present flow problems, it is necessary to include a 
large domain to get to the boundary conditions that are given at places 
far from the blades. The use of higher -degree polynomials will greatly 
reduce the niimerical work. 

If a number of cases are to be solved for a given set of blades 
on a given surface, it is advantageous to solve the problem on a large- 
scale digital computing machine. If a very high-speed digital machine 
is available, the simultaneous equation is best solved by Llebmann’s 
iterative process, which is the most simple to set up (reference 24) . 

If only a relatively slow-speed machine is available, the matrix 
process suggested in reference 18 is most suitable. 
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Matrix solution . - In the matrix solution, the entire equation (8) 
Is used with p being kept constant during any one cycle. Several 
cycles are necessary for compressible flow because p Is In a continual 
process of change during each calculation according to equation (?). 

The matrix method Is found most satisfactory when several solutions are 
required for a given blading configuration, (for Instance, In the case 
where the flow Is to be determined through a blading at several inlet 
Mach numbers or inlet angles) because the factorization of the matrix 
is the most time-consuming part of the problem and once completed the 
matrix may be used again and again as one of the inlet flow conditions 
is varied. For example. In a problem similar to the cascade problem 
in which 400 grid points and fourth-degree polynomials were used. It 
took about 60 hours to factorize the matrix on an IBM Card Programmed 
Electronic Calculator, and solving the set of values of ijf at all 
grid points for a given set of J values only takes 2 hours. The 
rounding-off error is found to be very small. 

Relaxation solution . - In the relaxation solution, the entire equa- 
tion (8) is used in the calculation of the residuals; but only the first 
two terms are used in the actual relaxation. As in the matrix solution, 
several cycles are required for compressible flow because p is not 
allowed to vary during the relaxation process whereas actually p is 
changing according to equation (7). Consequently, after each cycle 
new values of p are calculated and used in the determination of an 
improved J as in the matrix method. 

Iteration method of Llebmann . - The iteration method of Liebmann 
is simplest to set up on a digital computing machine; but it is the 
slowest method, and therefore it can only be used for the present 
problem on a high-speed machine. The boimdary conditions and the 
between-cycles J correction process are the same as those in the 
matrix and relaxation methods. 


ACCURACY OF FINAL RESULTS 

Several factors affect the accuracy of the final results: the 

degree of Lagrangian polynomial used, the size relative to the airfoil 
of the network spacings, the smallness to which the residuals are 
reduced, and the number of cycles carried out for compressible flow. 

The fourth-degree Lagrangian polynomial has been found to be the 
most applicable to the type of calculation under consideration. The 
use of a second-degree polynomial, in general, requires a smaller 
grid spacing, necessitating a much larger number of network points than 
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the fourth-degree polynomial. Polynomials of higher degree than the 
fourth are not desirable within the channel, because fitting the points 
necessary for these higher-degree polynomials circumferentially inside 
the channel may require a smaller grid spacing. Also the accuracy near 
the leading and trailing edges of the blade is best obtained by using 
small spacing. Outside the channel, however, the use of a higher- 
degree polynomial would not only reduce the amount of work but also 
would specify more completely the condition that the flow repeats cir- 
cumferentially every pitch angle. 

Formulas that can be used, after the problem has been completed, 
to determine the order of magnitude of the error at any point are given 
in references 18 and 23. These formulas give only qualitative informa- 
tion concerning the obtainable accuracy for the grid size and the 
degree of polynomial chosen, especially when there is no solution of a 
similar problem available. Because no solution is available with which 
to apply these formulas at the start of the example to be discussed in 
the following section, a somewhat smaller grid spacing than that 
required for the use of a foiurth-degree polynomial across the channel 
is chosen. (See fig. 3.) Results obtained in that example indicate 
that the grid chosen gives sufficient accuracy for the present p-urpose. 
If a more acc\n*ate result is desired around the nose or tail region, 
either small spacing can be used in these regions or the flow in these 
regions can be improved by using a finer net after the solution is 
obtained with the original coarse net. 

For the matrix and iteration methods, the final acciiracy is 
determined by the choice of grid size and degree of polynomial; for 
the relaxation method, however, the degr-ee of residual reduction must 
also be included in an accuracy determination. The smallness to which 
the residuals can be reduced depends essentially on the amount of 
labor put into the problem; but the smallness may be limited by the 
use of too large grid spaces or too low a degree of polynomial. An 
approximate method for estimating, at the beginning of the problem, 
the degree of residual reduction necessary is outlined as follows: It 

is assumed that an error of 1 percent in addition to that introduced 
by the finite -difference equations will be tolerated in the first order 
derivatives of For example, the following eqmtion can be written 

in the f?-direction by the use of the usual differentiation formula: 



( 12 ) 
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where the ijf's are the correct values of ijr, the E's are the mag- 
nitudes of the errors in &cp is the grid spacing in the cp- 

direction, and the superscripts c, d, and e refer to three eq.mlly 
spaced grid points (fig. 3). Ecuation (12) may be rewritten as 

^ E° - E^ ( 13 ) 

^cp j ^ 26 (p 25 cp 

J.C _ -pe \lrc _ ,j<e 

where, in this case, — 55 is less than 1 percent of — 2QZ'~ 

specified. As an example, take equal to 1000, then 

E° - E® < 10 . Because a constant magnitude error is maintained through- 
out any one area in the network, |e®| = 1e®| <5. The final residual 
must be such that a change of \lf of no more than 5 is necessary to 
reduce the value of the residual to zero. The absolute value of the 
residual can therefore be no greater than the absolute value of 5 times 
the relaxation coefficient at the point in question when the final 
accuracy is reached. For relatively small grid spacings, this method 
should give a good estimate of the magnitude of allowable residual. 
Because the relaxation coefficient at the point in question varies from 
point to point, it is advisable to place the corresponding maximum 
allowable residuals on the relaxation sheet for ready reference. 


NUMERICAL EXAMPLE 

* 

The method described is applied to analyze the detailed flow past 
t\irbine blades on a cylindrical surface as shown in figure 3. In such 
a case, a equals 0 , I becomes z, and it is more convenient to 
use a distance on the circinnference y(=rcp) instead of the angular 
coordinate cp (the reduced equations are given in the appendix) . 

Because only one Incompressible and one compressible solution were to 
be obtained at the design condition, the finite-difference equations 
were solved by the relaxation technique. This particular cascade was 
chosen because of the relatively high turning and the considerable 
thickness of the blade involved, thus giving the method described a 
thorough test. Experimental data for several inlet Mach mmibers are 
available for this cascade, increasing its value as a problem choice. 

The pitch and axial chord of the cascade are 1.017 and 1.5 inches, 
respectively. 

The grid, shown superposed on the cascade in figure 3, was obtained 
by dividing the pitch into eight equal parts, thereby arriving at a grid 
spacing of 0.1271 inch, which is used in both y and z directions. 

As suggested, fourth-degree Lagrangian polynomials were used throughout 
the domain to obtain the finite-difference equations. 

The reference lines outside the channel were drawn parallel to the 
flow at the inlet and the exit. The inlet angle is given as 41°18' . 
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The exit angle, as computed hy the measured pressure far dovrastream 
and the continuity equation, is equal to -52°57 ' , which checks quite 
well with the inverse-sine empirical rule. The fixed-angle method was 
used to acco\ant for the inlet and exit hoimdary condition. The inlet 

and exit houndaries were placed a distance equal to 1^ times the pitch 
from the cascade channel, which was found to be more than sufficient. 

The results obtained from the relaxation solution of this cascade 
serve a threefold purpose: (l) They show the usefulness of the method 

discussed) (2) they give some detailed information concerning the var- 
iation of flow in a two-dimensional cascade that may serve as a basis 
for devising some simple, good approximate methods of analysis) and (3) 
they give some useful information for the through-flow analysis in 
turbines having thick blades. 

The results obtained in the incompressible solution are shown in 
figures 5 to 13. The variations of the magnitudes of W 2 , Wy, and W 
across the channel are shown in figures 5 to 7 . The variation in Wy 
is seen to be larger than that in and larger at the leading por- 

tion of the channel than at the trailing portion of the channel. The 
variation of the velocity components across the channel can be accu- 
rately represented by a second-degree function in the y-direction. 

Several streamlines, with equal mass flow between them, are shown 
in figure 8. Also shown is the mean channel line and the mean camber 
line. The mean streamline, which divides the mass flow in the channel, 
is seen to have a cttrvature less than either the mean channel line or 
the mean camber line and is closer to the suction surface than the 
pressure surface. (The mean camber line has the highest curvature.) 

The slope of the mean streamline is compared with those of the blade 
mean camber line and the mean channel line in figure 9. Inside the 
channel the absolute value of the mean streamline slope is found to 
be lower than those of the mean camber and mean channel lines. Pro- 
ceeding outward from the inlet channel proper the mean streamline 
slope at first increases above and then approaches the specified inlet 
and exit values. Figure 9 shows that the specified values are reached 
approximately one pitch distance from the blade . 

The variations of the mean streamline velocity components and their 
derivatives with respect to z are shown in figures 10 to 13. In fig- 
ure 10, the variation of W 2 along the mean streamline is compared with 
the variation in channel width. Inside the channel, the increase in 
Wz due to reduction in channel width is on the average about 4 percent 
higher than that given by the one -dimensional calculation based on the 
reduction of the channel width. The difference is due to the nonlinear 
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variation of W 2 across the channel. The influence of the blade thick- 
ness also extends a short distance both upstream and downstream of the 
channel. With a finite blade thickness and a finite blade spacing^ the 
tangential velocity on the mean streamline (fig. 12) is seen to rise 
above its inlet value in front of the leading edge, to decrease rapidly 
for the first half of the channel, and then to decrease slowly to a 
value below its exit value for a short distance downstream of the trail- 
ing edge. The derivative with respect to z of the tangential velocity 
is shown in figure 13. 

The compressible solution is obtained for an inlet Mach nvunber of 
0.42 in order to compare it with the available experimental data at 
that inlet Mach number. The results obtained (figs. 14 to 23) are pre- 
sented in the same manner as those for the incompressible case. The 
velocities (figs. 14 to 16) show a somewhat higher rate of variation 
across the channel and can again be approximated by second-degree func- 
tions. 

The difference in streamlines between the conipressible (fig. 17) 
and the incompressible (fig. 8) solutions is appreciable. In general, 
the streamlines are pushed farther away from the suction surface. As a 
result, the mean streamline now approximates the position of the mean 
channel line better than in the incompressible case but still has a 
significantly lower curvature. 

The comparison between the slope of the mean streamline and those 
of the mean camber line and the mean channel line in the compressible 
case (fig. 18) is similar to that in the incompressible case (fig. 9) . 
Instead of the axial velocity in the incompressible case, the specific 
mass flow pW^ on the mean streamline is compared with the channel 
width in the pitch direction in figure 19. The comparison is similar 
although the difference is somewhat greater in the compressible case. 

The tangential velocity on the mean streamline obtained in the com- 
pressible solution (fig. 20) and its derivative (fig. 21) are again 
s imi lar to the corresponding cinrves obtained in the incompressible 
solution. 

The velocity obtained around the blade in the compressible solu- 
tion is compared in figure 22 with the experimental data obtained at 
the NACA Lewis laboratory. The agreement is better at the pressixce 
swface than at the suction surface, as might be expected, and is sat- 
isfactory as a whole. Constant Mach number contours obtained in the 
relaxation solution are shown in figure 23. With an inlet Mach number 
of 0.42 and exit Mach nimiber of 0.55, the variation of Mach number from 
the pressure surface to the suction surface is quite large. The maximum 
Mach number on the suction surface exceeds 0.8. 
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SUMMARY OF RESULTS 

A method of obtaining the flow of a nonviscous compressible fluid 
past arbitrary compressor or turbine blades between two siirfaces of 
revolution is presented. The equations of continuity and motion 
obtained for such flow are combined into a nonlinear second-order par- 
tial differential equation in terms of a stream function defined for 
such flow. Numerical solution by the use of differentiation coeffi- 
cients for unequally spaced grid points is suggested. Means for satis- 
fying the boimdary conditions outside the channel and solutions by the 
relaxation method with manual computation and by a matrix method with 
a large-scale digital computing machine are described. 

Satisfactory results were obtained by the method described in the 
investigation of the detailed flow variation of a compressible fluid past 
typical high-solidity highly cambered thick turbine blades on a cylin- 
drical surface. The variations in fluid properties across the channel 
appeared to be representable by a second-degree function. The mean 
streamline approximately followed the shape of the mean channel line 
of the cascade and had lower curvatxare. The variation of specific mass 
flow along the mean streamline followed the trend of the variation in 
the channel width. In general, the variation in specific mass flow was 
significantly higher than that given by the ratio in channel width and 
the effect extended outside the channel. The velocity distribution 
around the blade obtained in the theoretical calculation compares very 
well with experimental values. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, February 5, 1951. 
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APPENDIX - DERIVATION OF EQUATIONS 

If the fluid flow in turbomachines can be assumed to occvir on sur- 
faces of revolution, it may be analyzed on a two-dimensional basis with 
the flow passing through a number of stream filaments of revolution 
(fig. 1(a)), the thickness of which may be taken as that obtained in a 
throu^-flow calculation (reference 22) . Such flow can be most con- 
veniently described by a set of orthogonal coordinates I and cp , 
where 7 is the arc length of the generating line of the mean surface 
of revolution in the meridional plane, and cp is the angle (fig. 1(b)). 
When the flow across an element defined by dZ and dcp is considered 
(fig. 1(c)), the continuity relation for steady relative flow gives 

pWj Trdcp + pWcp + I dZ 


dZ^ dcp - 

Factoring out dcp dZ and allowing dcp and dZ to approach 
zero in the factor yield 


pWz + 


B(pW^) 

dZ 



+ sin o 



§1 

dZ 



- pWjr 


dT ^(pW^) 

dZ ■ ~^§Z 


pWjTsin a - T 


5(pW(p) 

bcp 


= 0 


By the use of the relation 


dr 

dZ 


sin a 


equation (Al) can be written as 

3(TpW^r) a(Tp%) 


(Al) 


(A2) 


( 1 ) 


For steady adiabatic frictionless flow with uniform inlet condi- 
tions, the equation of motion in the circumferential direction can be 
obtained from equation (l4a) of reference 22 as 
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r 


1 



0 


or 


2 SWj SWqj Wqj sin a 


2ix> sin a = 0 


( 2 ) 


r Sep d1 r 


Special Case of Conical Siorface 


For "the special case of ‘two-dimensional flow be^tween "two neighbor- 
ing surfaces of revolution, where the mean stream surface can be taken 
as the surface of a right circular cone whose axis coincides with the 
z-axis, 0 becomes a constant. If I is measured from the apex of 
the cone. 


Alternative forms of equations (l) and the preceding two equations can 
be obtained in terms of 2 and an angle measured on the conical sur- 
face, as given in reference 17. 


r = Z sin a 


and equations (2) and (s) become, respectively. 


r Sep SZ 


1 SW^ Wbp 

r Sep SZ Z 


2 c 0 sin 0 = 0 


(A4) 


and 




(A5) 
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Special Case of Cylindrical Siarface 

For the special case of flow on a stream filament of revolution 
whose mean surface is a cylindrical sin:face, a becomes zero, X becomes 
z, and equations (l), (2), (4), and (s) reduce to the following forms: 


S(TpW^) S(TpWy) 

5z ^ 5y ~ ^ 

^z Sy ^ 



TpWy 




J 


(A6) 

(A7) 

(A8) 


S^ilr S In T St , /SlnpSf. Slnp St\ 

^ - -ST- SJ * ^ - V~5T- Si + Sr~ 


(A9) 


where y = rcp and r is a constant. When r is constant^ the rela- 
tion “between density and derivatives is simplified to 



(AlO) 


When the general table or graph constructed according 


is used, Z becomes simply 


\2 




and $ becomes 


to equation 


(7a) 



CS T2pj^i2) 


-1 
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(a) Stream filament of 
revolution. 


(b) Mean stream surface 
of revolution. 



(c) Element. 
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(d) Applications to axial-flow, mixed-flow, and radial-flow turbomachines. 

Figure 1. - Plow on arbitrary surface of revolution and stream filament of revolution. 
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Figure 2. - Values of J and k for two grid points near boundary using five-point system. 
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Pitch angle 


ro 


r 



Figure 3. - Reference lines upstream and dovmstream of blades and grid points. 


I9T2 


NACA TN 2407 


2151 


4 


MCA TN 2407 


25 


Inlet 




(b) Streamline-adjustment method. 


Figure 4. - Treatment of inlet boundary condition. 
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Figure 5. 


- Variation of axial velocity across channel for incompressible solution. 
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Figure 9. - Comparison of slopes of mean camber line, mean channel line, and mean streamline obtained in 

incompressible solution. 
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Figure 10. - Comparison between axial velocity along mean streamline and channel - 

width ratio for Incompressible solution. 
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Figure 11. - Variation of total derivative of axial velocity on mean streamline 
with respect to axis for Incompressible solution. 
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Figure 12. - Variation of tangential velocity on mean streamline 
with respect to axis for incompressible solution. 



Figure 13. - Variation of total derivative of tangential velocity on 
mean streamline with respect to axis for incompressible solution. 
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Figure 14. 


Variation of axial velocity across 
Inlet Mach number. 


channel for 
0.42. 
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solution. 
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FlKure 15. - Variation of magnitude of tangential velocity across channel for compressible 

solution. Inlet Mach number, 0.42. 
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Figure 16. 


- Variation of 


magnitude of resultant velocity across 
solution. Inlet Mach number, 0.42. 
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for compressible 
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Figure 18. - Comparison of slopes of mean camber line, mean channel line, and mean streamline 
obtained in compressible solution. Inlet Mach number, 0.42. 
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Figure 19, 
line in 


- Comparison between variation of specific 
compressible flow and channel-width ratio. 


mass flow along mean stream- 
inlet Mach number, 0.42. 
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Figure 20. - Variation of tangential velocity on mean streamline with 
respect to axis for compressible solution. Inlet Mach number, 0.42. 



Inlet Mach number, 0.42. 
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Figure 22. - Comparison of velocity distribution on blade obtained In 
numerical solution with experimental data. Inlet Mach number, 0.42. 
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Figure 23. - Mach number contour as obtained in the 
relaxation solution. Inlet Mach number, 0.42; exit 
Mach number, 0.55. 
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